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The development of Quantum Chaos in finite interacting 
Fermi systems is considered. At sufficiently high excitation 
energy the direct two-particle interaction may mix into an 
eigen-state the exponentially large number of simple Slater- 
determinant states. Nevertheless, the transition from Poisson 
to Wigner-Dyson statistics of energy levels is governed by the 
effective high order interaction between states very distant in 
the Fock space. The concrete form of the transition depends 
on the way one chooses to work out the problem of facto- 
rial divergency of the number of Feynman diagrams. In the 
proposed scheme the change of statistics has a form of nar- 
row phase transition and may happen even below the direct 
interaction threshold. 

PACS numbers: 05.45. -fb, 73.23.-b, 71.70.-d 



I. INTRODUCTION 

The investigation of chaotic properties of finite in- 
teracting quantum systems has now a few decades his- 
tory pl-M . The Random Matrix Theory, originally intro- 
duced for modelling this problem, adequately describes 
the statistics of complicated high energy excitations (see 
Ref. S for an extensive review). However, the concrete 
mechanism of developing the chaotic behavior in the fi- 
nite Fermi system with increasing of excitation energy 
still requires the further understanding. The interest in 
this problem has even increased in last years in partic- 
ular due to the recent explosion of activity in the field 
of Mesoscopic Physics (quantum dots, atomic clusters, 
helium droplets etc.). The experimental observation of 
electronic excitations in quantum dots |9[| has stimulated 
the progress jl^-Q in theoretical investigation of chaotic 
disintegration of single electron excitations in this sys- 
tems. The single particle excitations, being most easily 
produced experimentally, constitute an important part 
of Hilbert space of complicated system. However, the 
number of such states is exponentially small compared 
to the total number of excited states. The onset of chaos 
in the majority of many-particle states is the subject of 
this paper. 

We will consider finite, but sufficiently large Fermi sys- 
tems. The particles are supposed to occupy almost all 
orbitals below the Fermi energy ep of mean field single 
particle Hamiltonian |l5[. The averaged single particle 
level spacing is A <C £f- We will denote as simple states 
the single Slater determinant states having definite occu- 
pation numbers (0 or 1) of all orbitals. The simple states 



are mixed by the two-particle interaction and the ma- 
trix elements V are supposed to be small and Gaussian 
random 
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Here g may be interpreted as dimensionless conductance 
in the case of quantum dot |ll|Jl^ . This simple model 
was in fact introduced long ago in Refs. y,^. 

Due to the weak interaction Eq. (|l|), the low energy 
excited states are almost unperturbed. With the increase 
of energy £ the density of states grows rapidly and the 
exact wave functions at £ ^ A are formed by many 
chaotically mixed simple states. However, this chaotiza- 
tion of excited states proceeds in a very nonuniform and 
complicated way. For illustration we have shown on the 
Fig. 1 different energy scales relevant to this problem. 

The complexity of excited state is usually characterized 
by the inverse participation ratio (IPR) / = X]fc IV'iC"^)!^! 
where ipi{k) is the amplitude of fc-th simple excitation in 
the i-th exact wave function. Physically / is the effective 
inverse number of simple states mixed into exact one. 

While the density of single particle states vi = 1/A 
does not depend on energy, the total density grows ex- 
tremely fast with £ being I'totai ^ A~^ exp(27r-y/£/6A) 
|l6| . We use large £ for the total energy of the many- 
particle configuration and small e for the energy of single- 
particle orbitals. Due to the huge number of states (ana- 
log of the micro-canonical ensemble in thermodynamics) 
the averaged occupation number for given orbital for ar- 
bitrarily chosen "typical" simple state is given by the 
usual Fermi-Dirac distribution 



f{e - £f) = n{e) = exp 
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where T = V6f A/tt [0. Note, that we still have not 
taken into account the interaction V , but only explore 
the wide statistics of single-Slater-determinant states. 

The low energy thresholds on the Fig. 1 correspond 
to the decay of single quasi-particle. Quasi-particles in 
quantum dot may be seen as peaks broadened by the in- 
teraction (bunches of (5-peaks) in the one particle spectral 
density up to the energy e ~ Ag |10|. Above e ~ ^^/g 
the broadening arise due to a decay of quasi-particle into 
two particles and one hole and may be explained by the 
Fermi golden rule. Below this energy the width predicted 
by golden rule turns out to be smaller than three-particle 
level spacing. However, as it was shown recently p2[ up 
to energy e ~ A-\/g/ Ing the quasi-particle disintegra- 
tion still proceeds through the effective interaction with 



the many excited particle and hole states. A nonuniform 
distribution of the matrix elements of effective interac- 
tion in this case leads to a very unusual distribution of 
spacings between peaks and widths of the bunches |1J] . 
The original intention of this research was to apply 
the methods developed for description of quasi-particle 
disintegration for the investigation of the full spectrum 
of complicated system. However, the onset of chaos in 
the many-particle excitations provided us with a number 
of surprising new effects not present in the single quasi- 
particle case. Only the estimate of the effective high or- 
der matrix element in the section 3 may be considered as 
straightforward generalization of the analogous estimate 
in the Ref. |1J] . In the section 2 we will show, how the di- 
rect interaction between particles may mix into one wave 
function an exponentially large number of simple states 
even inside the region of Poissonian spectrum. In the sec- 
tion 3 we consider the self-consistent scheme of onset of 
chaos based on the analysis of the scaling behavior of the 
direct and effective interaction. However, in the following 
section 4 we propose an another scenario of the transi- 
tion to Wigner-Dyson caused by the factorial divergency 
of the high order interaction. The detailed comparison 
of the features of these two scenarios is given at the end 
of paper (Sec. 5). 



II. MIXING BY THE DIRECT INTERACTION 

For weak interaction the strong mixing of different 
quantum states may take place only if their energies oc- 
casionally turn out to be very close. Therefore, consider 
first only two states \i) and |/) with energy difference e 
connected by matrix element V. The simple calculation 
|p8| gives the IPR for exact wave function a\i) + (3\f} 



/ = a^ + r = 1 - 



2V' 



4^2 



(3) 



Let us fix the state |i) and consider a set of states |/) with 
a density dn/de — A7 . It is convenient to introduce the 
intermediate scale Eq; \V\ -^ eg <^ A/. With the proba- 
bility 1 — 2eo/Af there are no final states falling into the 
strip — £o < e < £o and / = 1 up to negligible correction 
~ V'^/AfEQ. With the small probability 2£o/A/ one of 
the states |/) enters the strip. In order to find / in this 
case one has to average the Eq. (^) over the interval of 
energy 
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and then average the \V\ over the ensemble of matrix 
elements. Combining the two contributions one finds 
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This simple formula is basic for our perturbative esti- 
mates. The averaging of the modulus of V for strongly 
non-Gaussian high order effective matrix elements will 
be the main subject of the next section. 

The goal of this paper is to consider the mixing of com- 
plicated excited states. Occupation numbers for these 
states are described by the Eq. (H). The density of di- 
rectly accessible states, or the density of pairs of occupied 
1, 2 and unoccupied 3, 4 states with £1+62 = £3+^4 which 
may be connected directly by the interaction is 
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Here /» = f{xi) (g), a;i,2 = £1,2 - £f, 2^3,4 ^ £f - £3,4 
and we have used that 1 — f{x) — f{—x). The factor 
(2!)~2 accounts for the identical initial and final parti- 
cles. It is easy now to use the Eq. (||) in order to find 
the IPR. However, this result may be even further im- 
proved. Eq. (H) with Af = A2 describes the situation 
when only for one two- by- two transition 1,2 — > 3,4 the 
energy difference turns out to be of the same order of 
magnitude with the matrix element. In the second or- 
der in V the leading contribution to / is given by double 
event 1, 2 ^ 3, 4 and 5,6^7, 8. However, as far as none 
of the orbitals 5 — 8 not coincides with any of 1 — 4 these 
two corrections are completely independent and corre- 
sponding contributions to / should be simply multiplied. 
Effectively, all one-particle states are chosen from ^ T/A 
orbitals around Ep- Therefore, if the number of couples 
of pairs falling into the energy interval '^ \V\ is small 
compared to ~ T/A one may neglect the interference 
between different direct transitions. Again we may intro- 
duce the intermediate scale £0; \V\ <C £0. The number of 
directly accessible final states is 2£o/A2 ^ 1, but until 
£o/A2 <C T/A the IPR is found simply as a product over 
these channels 



I = \ l-2n 



2£o 



= exp 



exp 




(7) 



Here we have averaged \V\ for the Gaussian random in- 
teraction Eq. @j. This result is vafid for A/g < A2T/A, 
or T <C Ay/g. The exponentiation of the IPR is a direct 
consequence of the two-particle character of the interac- 
tion. It would be very inconvenient to investigate such 
effects if one will consider the Hamiltonian of the system 
as a random (sparsed) Matrix in the full Hilbert space. 
Surprisingly, to the best of my knowledge, the simple 
effect Eq. (^ has not been considered yet e.g. for de- 
scription of the excitations in complex nuclei. 

Due to the Eq. (j^) at £ > Ag2/3 the number of com- 
ponents in typical wave function becomes exponentially 



large. However, within the range of vahdity of the Eq. (|^) 
only a small (also exponentially) part of close energy lev- 
els is mixed in any exact eigenstate. Therefore this simple 
effect is not enough for creation of full Quantum Chaos. 



lengthy discussion may be summarized in the formula 
for the probability distribution for the effective matrix 
element (||) 

^*'n^)-^^; 5-'<l^l/A<5-^ (10) 



III. HIGH ORDER EFFECTIVE INTERACTION The averaged modulus of the matrix element in this case 



Corrections to the Eq. (R) arise after taking into ac- 
count of the high order effective interaction between 
three-particle states, then four-particle etc. The example 
of corresponding tree- type Feynman diagram is shown on 
Fig. 2. The tree- type diagrams are maximally enhanced 
due to a wide statistics of initial and final orbitals. Con- 
tributions from the diagrams with closed loops are small 
like some power of A/T. Generalization of the Eq. (g) 
for the case of n-particle interaction gives 



1 

A^ 



HJ2^i) TT f ^ 
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i<2ri 



A 



i(2n)!A2« 
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Again (n!)^^ accounts for the identical fermions. Note 
that (2n)! ~ 2'^"{nl)'^ and therefore the integration over 
the energies Xi does not lead here to any n! suppression. 
This is in strict contrast with what happens for decay 
of quasi-particle ||lj] , where the main suppression of the 
density of final states i^n arose due to the integration over 
the energies of final particles. 

Sufficiently more complicated turns out to be the cal- 
culation of averaged modulus of the effective matrix el- 
ement entering both the Eq. (^) and the Eq. (|^). The 
simplest variant of the effective interaction is given by 
the three-particle interaction arising in the second order 
of perturbation theory 
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Formally this interaction mixes together three states |1), 
1 2), 1 3). The use of the concept of effective interaction re- 
quires that the admixture of the intermediate state |2) in 
the total wave function is small compared to those of |1) 
and |3). This means that V12, V23 ^ |ei — £2! ~ |£2 — eal- 
The contribution to averaged IPR from strongly mixed 
(almost degenerate) |1), |2) and |3) should be consid- 
ered separately, but this correction is small like l/\ng 
compared to the total contribution of the effective inter- 
action. On the other hand, we will be interested in the 
tail of distribution of the effective matrix elements. For 
Gaussian distributed V (|l|) the probability to find a large 
individual matrix element of the original two-particle in- 
teraction is exponentially small. The large matrix ele- 
ments appears if one term in the sum in Eq. (0) occasion- 
ally has anomalously small denominator |ei — £2! ^ A. 
All the other (random) contributions in the sum in this 
case may be neglected for In 5 :^ 1. The result of this 



is IK 



(3), 



e//l 



Aln.g/52 



Higher order effective matrix elements are obtained 
by simple adding of more intermediate states into the 
Eq. (P). The extensive discussions of the problems aris- 
ing while calculating the averaged value of this matrix 
element may be found in Ref. ||lj|. For example, the 
averaged modulus of the four-particle effective matrix el- 
ement for given Feynman diagram is given by 



V 



(4) 
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V12V23V3 
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2,3 



£2e3 



(11) 
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1 de2 d£3 

|£2£3|2A2A 



Here compared to the Eq. (0) we have made the change 
of variables £1 — £2 -^ £2 and £1 — £3 ^ £3. The pre- 
integral factor ~ (A/g)^ in the Eq. (O) arose after aver- 
aging of Gaussian distributed (|l|) direct matrix elements 
Vij I . The upper bound for the both remaining integrals 



in the Eq. (Till) is given by the single-particle level spac- 
ing |£2,3| < A. Due to the logarithmic character of the 
integrals the order of magnitude estimate of the bounds 
is enough. More interesting is the origin of the bounds 
at small £2,3. Analogously to what we have discussed for 
the three-particle interaction the formula (jO) describes 
the mixing of four states |1), |2), |3), |4) in the third order 
of perturbation theory. The interaction of the four states 
may be replaced by the effective interaction of only two 
|1) and |4) if the admixture of intermediate states |2) and 
1 3) is small. This requirement for example for state |2) 
means that |£2| ^ IV12I ~ A/g (weak mixing of |2) and 
|1)) and |£2| > |V23V'34/e3| ~ (A/5)V|£3| (weak mixing 
of 1 2) and |4)). Another two inequalities for the vector 
|3) are |£3| » A/g and |£3| > {A/g)'^ /\e2\- Combining 
all this inequalities one finds for the range of integration 
in the Eq. (K^ 



A/5< |£2,3| < A 



(12) 



These bounds allows one to find the averaged matrix el- 
ement with the ~ l/lny accuracy: |V"gf|| ~ A{\ng)^/g^. 
In the similar way one may consider the five-particle, 
six-particle, etc. effective interaction. Unfortunately, the 
system of inequalities describing the range of integratin 
for the intermediate energies turns out to be sufficiently 
more complicated in these cases. 

The formulas like Eq. (0) may be used to find again 
the distribution of the effective matrix elements. In par- 
ticular for the four-particle interaction 
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Here the range of variation of £2 Sind £3 is the same as in 



the Eq. (11) Found in this way the probabihty distribu- 



tion in the general case has the form 



p(")(y) 



A 1 



X 



(n) 



5-"+^<|y|/A<g-i, (14) 



where x^"'\V) is a piecewise polynomial function of the 
order n — 3 of two logarithms ln(|y|/A) and Ing. Again 
in the simplest non-trivial example of n = 4 one finds 



(4)_ rindngVA) , g-3 < |T/|/A < .9- 



X' ' = 



ln(|F|g/A) , 5-2<|F|/A<ff-i 



.(") 



(15) 



Also notice the huge range of variation of V^r i (14^ for 
large n. Averaging of the matrix element over the distri- 
bution Eq. (Oh leads to 



'eff\ 



|y;;il~ A(ln5)"-V5" 



(16) 



Due to a complicated domain of integration over the in- 
termediate energies it is difficult even to estimate the 
value of the overall numerical constant here. However, 
the important result of Ref. JlJ] is that at large n this un- 
known overall constant is neither large like n! nor small 
like 1/n!. The appearance of logarithms in the coeffi- 
cients of wave function was also observed in |ig|] . 

Like in Eq. ([7|) the corrections generated by the many 
particle effective interaction are naturally exponentiated 
at least for n < T/A. Thus combining A„ (|) with V^^j 
( [l6| ) we obtain instead of the Eq. (0) 



j 4%/3/f 

F{x) ^l + Y^Fn 




(17) 



In the case of short range screened two-particle interac- 
tion V{fi,f2) ~ 6{fi — r2) the first coefficient here is 
Fi = 24-\/2/7r/5. This result for I may be compared 
with the perturbative result for the IPR for disintegra- 
tion of the single quasi-particle in the isolated Quantum 
dot found in Refs. |l|-0 



/„ 



1 
\ng 



Hy) 



E'^"?/" 
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where y = (e / A)'^ In g / g and e is the energy of quasi- 
particle. By analogy with the usual IPR, the Iq is a 



sum of fourth powers of the amplitudes to find given k- 
th quasi-particle in the i-th exact state Iq = J^i l''/'i(^)l^- 
The next section will be devoted to the more detailed in- 
vestigation of the possible properties of the scaling func- 
tion F{x). For single quasi-particle the coefficients 0„ 
decrease like nl~^ for large n, which makes very improb- 
able the disintegration of quasi-particle through dclocal- 
ization transition in the Fock space. 

The simple counting of the powers of £/A and g in Eq. 
(O) leads to the following scenario of the development 
of chaos: 

a). For A52/3 <£ < Ag/\ng the Eq. (0) is valid. The 
number of components in typical wave function is expo- 
nentially large, but described by simple multiplication of 
independent direct transitions. 

b). At £ ^ Ag/ Ing all terms in the series Eq. ( [l7| ) en- 
ters the game. The wave function is formed by a fractal 
combination of states distant in the Fock space [gO| . If 
the resumed scaling function F will be a regular function 
of the argument x (like it is for the cj){y) ([l8|)), the num- 
ber of components of the wave function will be small com- 
pared to the total density of states and the level statistics 
will be close to Poissonian. 

c). Starting from £ ^ Ag all quasi-particles constitut- 
ing the wave function have enough energy to decay via 
usual Golden-rule. Only here the Wigner-Dyson statis- 
tics is evident. The more formal criteria for onset of 
chaos will be the requirement that the product lEvtotai 
became of the order of ~ 1, where vtotai is the total den- 
sity of states described above the Eq. (pi). This condition 
is equivalent to the equation for unknown function F{x) 



xF{x) 



Q^/2 



In 5 



(19) 



Solution of this equation gives the energy of transition 
to rigid spectrum somewhere within 5/lng < E/A < g. 
Due to the exponential character of / and vtotai most 
likely this will be the narrow phase transition. 



IV. THE iV! SCENARIO 

Thus the power counting for the perturbative IPR 
Eq. (nrh may provide us with the formally self-consistent 
scenario of developing of chaos in finite Fermi systems. 
Now we are going to show, how taking into account of 
the high order behavior of the coefficients Fn may lead 
to the important revision of this straightforward scheme. 
Unfortunately, we are not able to perform the complete 
estimation of the asymptotics of F„. Nevertheless, the 
catastrophic change of the mechanism of onset of chaos, 
which we will consider, may partly compensate the lack of 
rigor. Therefore, even though based on certain assump- 
tions, the result of this section should be considered as 
one of the main results of the paper. 



The calculation of any term of the series in Eq. ([11 
consists of two main steps. First is the calculation of 
the density of final states, which we were able to per- 
form explicitly (g). The second step is the estimation 
of the effective matrix element. The matrix element for 
a given Feynman diagram is described by the Eq. (16). 
As we have told, this matrix element could not contain 
any n\ or n\~^. However, for many-body problem and 
high order of perturbation theory the same initial and 
final configurations may be connected by many different 
Feynman diagrams. We have shown on the Fig. 2 the 
example of corresponding tree-type n-particle diagram. 
The total number of diagrams like on Fig. 2 is 



^diagrams = Nd = n"' [n — l)!n! 



(20) 



Here n\ is the number of permutations of 1' — n' final par- 
ticles, (n— 1)! is the number of replacements of (screened) 
Coulomb interaction lines and n"~^ is the number of 
[n — f )-segment trees connecting n points 1 — n |21 1. Eq. 
( P0| ) gives us the number of diagrams of Schrodinger per- 
turbation theory. Among the Nd trees there are many 
diagrams having the same set (up to permutations) of 
the matrix elements of two-particle interaction Vij and 
different energy denominators. As it was pointed out in 
Ref. Q , these diagrams could not be considered as sta- 
tistically independent, which results in lowering of the 
mixing by the high order effective interaction. The aver- 
aged modulus of the effective matrix element found for 
single Feynman diagram (n6f) should be modified in the 
case of many interfering diagrams. However, due to a 
huge range of variation of the Ve// in the power law dis- 
tribution Eq. (Ill) the average value is saturated by a rare 
very large fluctuations of the matrix element. Therefore 
for n not too large the generalization of the Eq. ( [l6|) for 
many diagram case is achieved by the simple multiplica- 
tion by the number of diagrams (we will return later to 
the discussion of very large n). Finally, combining the 
Eqs. d), (|l|) and ^ we find 



Fn<n\ 



(21) 



We see, that the series in Eq. (|l7| ) most likely is the 
asymptotic series. Due to the unsatisfactory estimate of 
the number of diagrams and even the value of individual 
diagram we were able to draw only the upper bound for 
the Fn- However, it seems very unlikely that the cancel- 
lation between the diagrams of the Schrodinger pertur- 
bation theory will remove the factorial divergency of the 
coefficients of the series. In the following (if the oppo- 
site not indicated explicitly) we suppose that Fn ^ n\. 
The asymptotic series, common in the Field theory and 
High Energy physics, usually do not cause serious trou- 
bles (see for review Ref. p2i). In all known cases one 
simply breaks the summation on the smallest term of 
the series and use this smallest term also as the order of 
magnitude estimate of the rest nonperturbative part of 



the sum. The more refined strategy will be to perform 
the Borel summation, but this procedure (up to the same 
nonperturbative corrections) is equivalent to the break- 
ing of the series at the smallest term. In our example 
this means that one should break the summation in the 
Eq. ( p^ ) at no ~ g/S./h\g£ and the non-perturbative 
ambiguity in F should be 



<5F„ 



exp 



\ng £ 



(22) 



with some a r^ 1. However, in the known field-theoretical 
examples there exist the clear physical mechanisms (In- 
stanton, Renormalon) which allows to understand the 
origin of non-perturbative corrections and the necessity 
of breaking the series at the smallest term. In our many- 
body quantum problem we could find no physical moti- 
vation for the appearance of nonperturbative corrections 
of the kind (ra) and we see no reasons for breaking the 
summation in the Eq. (UT) at the smallest term. Of 
course, the Eq. ( p7[ ) is not exact, but corrections to 
F{x) are small like (Ing)^^ or A/£. Taking into account 
such corrections will help us, to sum up the original di- 
vergent series. Nevertheless, in order to make sense of 
our perturbative analysis of the onset of chaos in finite 
Fermi systems we should find the physically motivated 
cutoff for the n! divergent sum. The natural cutoff for 
Fn, which we see, is given by the total number of par- 
ticles which may enter into the diagram of effective in- 
teraction Uyjiax ~ T/A ~ ^j£ / A.. If Umax > x~^ the 
summation now goes above the smallest term of the se- 
ries. The consequences of the decision to sum up the se- 
ries above the smallest term are crucial. In the standard 
scheme the function F{x) is described by perturbation 
theory at small x and became completely nonperturba- 
tive (although it still may be smooth and finite) at a; '--^ 1. 
Within the new scenario F{x) blows up at 



-^n™„-.. ^ 



(n^a^lngf /gA)""°- - 1 , (23) 



which 

gives the chaotization threshold E^ ^ ^(s/ln^)^^'^. This 
estimate was done for the most strongly divergent series 
allowed by the Eq. (pi|). For example, if F„ ~ {n/2)\ the 
same calculation will lead to Ec ^ A(g/lng)'*/^. There- 
fore, the more accurate conclusion from the Eq. ( pi| ) is 
that due to the effective interaction including all ^ T/A 
excited particles the series for IPR Eq. (0) blows up at 
some (g/ In 5)2/3 <Ec<g/lng^. 

The two important peculiarities of the emergence of 
Quantum Chaos follow from this "factorial" scenario: 

I. The direct interaction connect only the states close 
in the Fock space pG]. Therefore, the Eq. (M) describes 
mixing of only a small part of close levels. On the other 
hand, the highest order effective interaction entering the 
Eq. (E3p mixes the huge number (the majority of all 
many-particle excitations) of very distant simple states. 



Thus, divergency of F{x) due to the n! leads to the tran- 
sition from Poisson to rigid (Wigner-Dyson) spectrum at 

£^ Se- 
ll. Due to the large power of x in the Eq. ( p3| ) the 
growth of F{x) with increase of energy takes place at a 
very narrow interval AS /£c ^ ^/nmax- This means that 
the change of statistics has a form of phase transition 
with relative width '^ l/umax- 

Our description of the development of chaos is suffi- 
ciently based on the estimate of Ti„iax (Ph) ■ ^^ '^^^ would 
be able to find another reasons for breaking the series 
at a term parametrically smaller than Umax ^ T/A the 
result for £c will be also different. For example, the new 
cutoff may arose due to a huge number of interfering dia- 
grams Eq. (po| ) . The joint probability distribution for the 
sum of many diagrams with power law individual distri- 
butions described by the Eq. ([l4| ) has a Lorcntzian form 
tejl with the width 7 having the form 



Slater determinant. The pure direct interaction in this 
case gives 



P(V) 



■• X(V).^^^,^,A». 



^3r! 



V2 + ^2 ^(^) 



(24) 



Here we have used for the estimate of the number of 
independent diagrams the upper bound Eq. (EG). This 
distribution is valid for \V\ < A/g. The function x was 
defined in the Eq. @. The factor xiV)/xi'l) allows to 
preserve the asymptotics Eq. ( p^ at \V\ '3> 7. Due to 
the logarithmic character of x this factor correctly ac- 
counts for the deviation from the pure power law in the 
distribution of single matrix elements. We see that af- 
ter taking into account of the interference of different 
diagrams the power law distribution of the total matrix 
element becomes strongly violated at \V\ ^ 7, not at 
\V\ - Ag-"+i as in the Eq. 0. The value of F„ found 
with the distribution P{V) (|24| ) becomes parametrically 
different from the estimate Eq. (pi| ) only at n'^^^ ~ g^^^ 
(even at nj„^j. ~ {g/ IngY^^ if one considers seriously the 
high powers of logarithms in the Eq. (^J)). However, this 
new variant of the cutoff leads to the same threshold for 



the onset of chaos as we have found before in Eq. (23). 
This coincidence of the results of two different approaches 
to estimation of the rimax rnay be considered as an indi- 
rect support for the transition to Wigner-Dyson statistics 
at£,^A(ff/ln5)2/3. 

At least some of the thresholds shown on the Fig. 1 may 
lie above £c- The physical meaning of these perturbative 
thresholds becomes less clear. Nevertheless, even above 
£c one may try to look for effects generated by the di- 
rect or low order effective interaction. The highest order 
effective interaction becomes important due to the huge 
combinatorics Eq. (20), but the corresponding matrix 
element is in general very small (up to 1^ ^ /^y^^Wmax^ 
Thus this high order interaction may be invisible (or at 
least suppressed) at small time in time-dependent prob- 
lems. For example, one may consider the time evolution 
of the wave packet ^(i), starting at t — from single 



l(*(0)|*(t))|^ 



~rt . 



r = 2V6Ag~^ {£/A) 



3/2 



(25) 



The formal power counting (like we have done below the 
Eq. (p^) gives for the range of validity of this result 
Ag^/^ ^ 5 <C Ag/hig. However, due to the logarith- 
mic behavior of all integrals (see e.g. the Eq. (pl|)) all 
values of the high order effective matrix elements are 
equally important within 5""+^ < \V\/A < g~^ (pl|). 
Therefore, the range of validity of the Eq. ( p5| ) may be 
changed strongly after taking into account of the contri- 
bution from large high order effective matrix elements. 



V. DISCUSSION 

The goal of this paper was to consider, how the chaotic 
mixing of nonintcracting eigenstates develops with in- 
creasing the excitation energy in finite weakly interact- 
ing Fermi systems. Measured in the units of the single- 
particle level spacing A the energy corresponding to the 
onset of rigid spectrum may depend on only one parame- 
ter g (|l|). Nevertheless, even with only one parameter we 
were able to introduce a number of different energy scales, 
shown on the Fig. 1 , which correspond to (or may corre- 
spond to) the different stages of development of chaos. 

The low lying thresholds on the figure e/A ~ 
(5/ In g)^/^ and e/A ~ g^^^ are associated with the decay 
of single quasi-particlc. These thresholds were considered 
in the Rcfs. pO|-p^ and their description was not the 
aim of this paper. However, the methods introduced in 
Refs. |l3,H for the investigation of below golden rule de- 
cay provides us with the useful tool for the investigation 
of onset of Chaos in complex many-particle excitations 
also. 

Technically our consideration is based on the two sim- 
ple equations (0) and M) for the fPR. The two non- 
interacting states are strongly mixed by the interaction if 
the corresponding energy difference turns out to be of the 
same order of magnitude with the matrix element. As a 
result, all the effects which we consider are determined 
by the product of the averaged modulus of the matrix 
element and the density of accessible states. 

The emergency of chaos in the majority of complicated 
states with many excited particles and holes was the main 
subject of this paper. This increase of chaos with the 
growth of excitation energy proceeds in a very compli- 
cated way, both technically and even more logically. In 
the simplest variant, mixing of the simple (single Slater 
determinant) states is caused by the direct interaction. 
The mixing by direct interaction opens at £" > Ag^^^ 
p5[ . Due to a large number of excited particles the in- 
teraction may proceed through a large number of inde- 
pendent channels. This results in the exponentiation of 



the IPR, as we have shown at the end of the Sec. 2. Thus 
even the direct interaction may mix together an exponen- 
tially large number of Slater-determinant configurations. 
However, this number is still small compared to the total 
number of levels of complicated system and the direct 
interaction could not cause the transition from Poisson 
to Wigner-Dyson statistics of energy levels. 

An alternative to direct interaction is the transition 
via effective three-particle, four-particle, etc. interaction 
considered in Sec. 3. The effect of this high-order in- 
teraction is much weaker than that of the direct inter- 
action at low energies, but they became comparable at 
£ ~ Ay/ In g. Above the energy £ ~ Ay/ Iny our pertur- 
bative approach is no more valid and we were forced to 
introduce the new unknown scaling function f (Ingg/gA) 
into the exponent describing the IPR in the Eq. (17). 
However, if one expects that F{x ~ 1) ~ 1, the logarithm 
of the number of components in the typical wave function 
aX £ r^ Ag/ In g will be in In g 3> 1 times smaller than 
the logarithm of the total number of states. This means 
that although the mechanism of mixing is changed and 
the number of components in the wave function became 
to grow much faster at 5 ~ Ag/ In g, this threshold does 
not correspond to the transition Poisson- Wigner-Dyson. 

The Wigner-Dyson statistics becomes evident only at 
£ ^ Ag. This energy is analogous to the golden rule 
threshold e ^ Ag^l"^ for the decay of single quasi-particle. 
Roughly speaking above the £ ~ Ag any particle from 
the strip ~ T around ep may decay via the golden rule. 

Formally, in order to understand, how the change of 
statistics proceeds between £ ^ Ag/ In g and £ ~ Ag 
one should analyze the high order behavior of the coeffi- 



cients Fn in the expansion of the function F{x) Eq. (17). 
Surprisingly, this analysis in the Sec. 4 instead of such 
understanding led us to the crucial revision of the con- 
sidered above straightforward scheme. 

First of all, as we have shown in the Sec. 4, the se- 
ries in the Eq. ( [l7| ) most likely is the asymptotic series. 
Although (|2l|) Fn ~ n\ is only the upper bound for the 
Fm it seems very unlikely that the factorial divergency 
of the coefficients may be washed out. Moreover, even if 
Fn ^ {vn)\ with any < z^ < 1, the results of the Sec. 4 
will be only slightly modified. The factorial divergency 
of the series would not be dangerous if one will sum it 
up e.g. via some variant of the Borel method. However, 
within the Borel prescription one is immediately faced 
with the problem of strange nonperturbative corrections 
like in the Eq. (p2h (moreover that Fn are positive by 
construction and the series is non-Borel summable). Be- 
ing unable to find the physical motivation for the Borel 
summation, we were looking for other way to make sense 
of the divergent series. The coefficient F„ is given by 
the sum of n -|- 1-st order tree type Feynman diagrams. 
Therefore, the natural cutoff for the series is given by 
the largest diagram which may be drown with available 
~ T / A particles. Namely nmax ~ T / A (one more way 



to determine Umax considered in the Sec. 4 leads to the 
same physical result). Within this scenario the rimax and 
the scaling parameter x = \ng£/gA in the Eq. ( p^ ) are 
two independent variables. At a; > l/rimax summation 
in the Eq. (O) goes above the smallest term. However, 
up to the energy £c ~ A(g/lng)^/^ all Umax terms of 
the series are small. At £ — £c even the broken series 
blows up and above £c our perturbative approach is no 
more valid. Because of the growth of F{x) at £ ~ £c 
proceeds mostly due to the last terms of the series or 
the most complicated diagrams of the highest allowed or- 
der of perturbation theory, this threshold corresponds to 
narrow phase transition from Poisson to Wigner-Dyson 
spectrum. We have attempted at the end of the previous 
section to find remnants of non-chaotic behavior above 
£c by considering the short time processes. However, this 
issue requires a separate and more careful investigation. 
To summarize, we have considered in this paper the 
two scenarios, how the quantum chaos may develops in 
the Finite Fermi systems depending on the way one will 
use for the summation of factorially divergent contribu- 
tions from the high order effective interaction. 

• If the series is summed up within the Borel pre- 
scription (i. e. broken effectively at the smallest 
term) we were able to find a few distinct stages 
of the chaotic behavior. Below the energy ^ Ag^/"^ 
the majority of eigenvectors is almost unperturbed. 
Between ~ Ag^/^ and ^ Ag/ In g direct interac- 
tion mixes together the exponentially large num- 
ber of simple states (M). However, the IPR is under 
theoretical control in this region and the statistics 
is clearly Poissonian. Above Ag/lng statistics of 
the energy levels is still Poissonian, but the IPR 
is described by some (unknown) scaling function 
F{x) (p7|). Transition to Wigner-Dyson takes place 
somewhere between £ '^ Ag/ In g and £ ~ Ag (see 
the Eq. (|l9|)). The transition may be relatively nar- 
row A£/£ ^ (Ing)"^ 

• In the second variant the cutoff of the series in the 
Eq. ( |l7| ) is determined by the number of the par- 
ticles in the largest allowed Feynman diagram and 
is independent of the scaling variable x. In this 
case the transition to Wigner-Dyson takes place 
at £c ~ A(g/ln(7)^/'^ and is narrow like A£/£ ~ 
A/T^{\ng/gf/\ 

To our current understanding, the second scenario is 
more physically motivated. However, our argumentation 
in this part is not completely rigorous and one still may 
find the reasons for summation a la Borel of the tree di- 
agrams. The important disadvantage (or advantage?) of 
the second scheme is that if it is correct, it will be proba- 
bly the first example of the phase transition caused by the 
factorial divergency of the series. Finally, although both 
of our scenarios are reach of new physical predictions, we 



were able to make only the order of magnitude estimates 
of the considered thresholds. Therefore, any experimen- 
tal confirmation (within the true experiment, or by nu- 
merical simulations) of the mechanisms of development 
of chaos considered in this paper are very desirable. 
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FIG. 1. The different energy scales corresponding to chao- 
tization of single- or many-particle excited states. 



FIG. 2. The tree type diagram contributing to high order 
effective matrix element. The screened Coulomb interaction 
is shown by dashed lines. 



